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Abstract: Use of inclusive Mt2 distribution for sparticle mass determination is discussed. 
We define new parameters M^uiin) and M™ 2 od (min), which are a kind of minimum of 
sub-systerm Mt2 values. Their endpoints are less affected by initial state radiations. We 
demonstrate that both rriq and rrig can be extracted from the endpoints of the distributions 
in the wide region of parameter space expected in CMSSM. We also present a comparison 
with MxGen distributions. 
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1. Introduction 

The quest of physics beyond the standard model is now reaching new phase with the 
LHC experiments [1,2] and various new dark matter (DM) searches [3,4]. In models with 
conserved parity for stability of a dark matter, new coloured particles may be produced 
in pairs at the LHC and each of them decays into particles involving the stable dark 
matter. This gives missing transverse momentum to the events, which would be useful 
in separating the signal from QCD backgrounds. Especially, in supersymmetric (SUSY) 
models squarks and gluino are produced and decay into lighter particles involving the 
lightest SUSY particle (LSP) which is usually identified as the dark matter. The signature 
is events with jets (+leptons) + missing transverse momentum. 

Because the dark matter cannot be detected in the LHC detectors, masses of SUSY par- 
ticles cannot be reconstructed as resonances. Various techniques to determine the masses of 
SUSY particles have been studied (See Ref. [5] for a recent review.). Important kinematical 
variables are peak position of M e g distribution, the endpoints of invariant mass distribu- 
tions of jets and leptons [6-15], the endpoint of Mt2 distributions [16-30]. Among those, 
most clean channel is probably the invariant mass distributions of jll channel although the 
branching ratio of this mode is rather small. 
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In the early stage of the LHC, it is more useful to study the distribution with jet + 
Ermiss channel [23,24]. The number of signal events is rather small in the discovery phase, 
therefore if we require multiple leptons in the final state, we would not be able to obtain 
an event sample large enough to do analysis. For this problem, quantities that are less 
sensitive to the detail of the decay pattern are useful. 

In previous papers [23,24], we have defined an inclusive M?2, which is calculated for 
any SUSY events at the LHC using hermisphere method. Hemisphere method is the way to 
define two sets of particles in the event, and each set (called hemisphere) contains particles 
coming from the same parent particle decay in the limit that parent sparticles are highly 
boosted into back-to-back [31,32]. We have shown an interesting correlation between the 
endpoint of the inclusive Mt2 distribution and mass of the heavier of squark and gluino. 
We have also demonstrated that a subsystem Mt2, Mj^sub), is useful to determine the 
gluino mass for m g < rriq case [24]. 

One of the important issues for mass reconstruction is existence of initial state radiation 
(ISR) [26,33-36]. When we produce heavy particles at the LHC, there would be a hard 
ISR jet which could have px about the same order of the mass of the produced particle. 
The ISR jets may be involved among the jets for the invariant mass or Mt2 calculations, 
then the distributions are smeared. In Ref. [26], it is shown that the effect of leading ISR 
can be removed by minimizing u Mt2 after removing a jet" for several choices of a removing 
jet. Smearing due to ISR is significantly reduced by this procedure for a limited example 
gg production process followed by g — > XiJJ- 

In this paper we define a M^mm) to reconstruct the endpoint of Mt2 distribution 
under the effect of ISR for general squark and gluino production processes. The definition 
is also consistent with the previously defined M^2(sub) to obtain m g when m g < rriq. We 
extend our studies so that we can cover rriq < m g case as well, so that m g is reconstructed 
by the appropriate choice of quantities. For this purpose, we define a M™ 2 od (min) for the 
events that two highest pr jets coming from two body decays with large mass differences 
between parent and daughter particles and there are two jets which are prominent over the 
other jets. As its nature, the M^ od (min) works for rriq < m g . We provide parton level and 
jet level studies and demonstrate that these quantities are actually important. 

This paper is organized as follows. In section 2, we introduce basic parameters such 
as Mt2 and Mj-2(sub), which has been defined in the previous papers by using hemispehre 
algorithm, and summarize previous analyses. We also define new parameters M^uiin), 
and M™ 2 od (min). These parameters are defined to reduce ISR effects but for different 
squark and gluino decay patterns. Section 3 is for parton level study of hemisphere analysis, 
which is efficient for inclusive analysis. Section 4 describes the parton level distribution 
of M T2 (min) and Mf? 2 od (min) with and without ISR. The needs of M^ 2 od (min) will be 
explained in this section. In section 5, we discuss jet level Mt2, M^mm) and M™ 2 od (min) 
distributions. The measured endpoints depend on the dominant decay pattern of sparticles. 
We found the decay pattern can be determined from simpler quantities such as px of the 
two highest px jets and hermisphere mass distributions. The inclusive Mt2 study can 
be applied in the early stage of the LHC. The expected distribution at y/S = 7TeV and 
J Cdt = 1 fb _1 is shown in section 6. Section 7 is devoted for discussions and conclusions. 
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Figure 1: Total SUSY production cross sections for the model points in Table 1 as a function of 
mo for fixed gaugino masses at 7 TeV (left) and 14 TeV (right). The CMSSM points are produced 
by ISASUSY [50], and the cross sections are estimated by HERWIG [47-49]. The ul masses and 
names of the model points are also shown in the figure. 



2. Inclusive Mt2 distribution, ISR, and mixed production 

The discovery of super symmetric particles at the CERN LHC has been studied in depth 
[37-46]. In recent articles [45,46], the discovery region at \fs = 7 TeV and J Cdt = lfb -1 
in CMSSM has been studied. One can read from Fig. 13 in Ref. [46] that the region 
tuq < 850 (750) GeV at M 1 1 2 < 210 (250) GeV can be explored, respectively. Here mo and 
Mi/2 are universal scalar and gauino masses. In Fig.[j], we show the total SUSY production 
cross section as a function of m for My 2 = 210 (250) GeV at y/S = 7 TeV (left) and 14 
TeV (right). The model points appeared in Fig.[l] are defined in Table [l[ The cross section 
at each model point is calculated by HERWIG [47-49]. The cross section at the boundary 
of the discovery region for y/S = 7 TeV with J Cdt = 1 fb _1 is roughly asusY = 1 pb. After 
the 7 TeV run, the collision at 14 TeV is scheduled. The cross sections at 14 TeV is about 
10 times larger than those at 7 TeV. 

By studying production cross sections of sparticles at the LHC, patterns for sparticle 
masses can be investigated. The production cross section is sensitive to the mass scale of 
produced particles, as can be seen in Fig. [l], and it is an important piece of information on 
the mass scale. Increasing mass by 100 GeV leads 50 % reduction of the cross section at 
rrig ~ rriu L ~ 500 GeV. The uncertainty of the squark and gluino production cross section 
would be less than 10% [51]. 

One can determine the mass spectrum directly from kinematical distributions of spar- 
ticle decay products. If both squark and gluino masses can be determined, one can further 
study their interaction by comparing the measured production cross section with the the- 
oretical value, because the dependence on the other parameters are rather small [52,53]. 

In SUSY models with conserved R parity, two LSPs escape from detection for each 
SUSY event, therefore masses of sparticles cannot be measured as resonances. (In this 
paper the LSP is assumed to be the lightest neutralino, Xv) Endpoints of invariant mass 
distributions of jets and leptons distribution can be used to determine the masses. For 
example in q — > x% ~ ^ I ~ * Xi produce a jll final state and endpoints of the invaiant mass 
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Table 1: Masses and important branching ratios in some CMSSM points. Points 1 — 5 for fixed 
gaugino mass Mi/ 2 — 250 GeV, Points 1' — 5' for Mi/ 2 = 210 GeV. Mass spectrum are calculated 
by ISASUSY [50] and cross sections are calculated by HERWIG [47-49]. 

distribution mju, mji(min), m,ji(max), m\\ can be used to determine the all sparticles 
masses involved in the decay. However, branching ratios of favorable modes such as jll 
are generally small, and they may not be useful in the early stage of the experiments. 
Therefore it is important to find inclusive quantities sensitive to squark and gluino masses 
so that one can use the most of the signal events at the LHC. 

Typical energy scale of sparticle production processes may be estimated from the 
"effective mass", 

M cS ^Y,\PTi\ + Ef s \ (2.1) 

i 

where in this paper the sum is taken for the jets with px > 50 GeV and r] <2.5. A peak 
value of the M c r distribution is correlated with squark and gluino masses, although the 
relation is rather qualitative [54]. 
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More recently, we have proposed an "inclusive Mt 2 v which is calculated from jet (and 
lepton) momenta and the missing transverse momentum [23,24]. The definition of Mt 2 is 
the following [16,17]: 

M T2 = min [max(M T (p Vl ,pX),M T (p V2 ,pX))}. (2.2) 

Here v\ and v 2 are two selected visible systems in an event, and pf is a test LSP momentum 
with a test mass m x . If v\ and v 2 arise from two on-shell particles A\ and A 2 as Ai — > ViXi, 
respectively, the Mt 2 is bounded from above as 

Mt 2 < max(m J 4 1 , mA 2 ) for m x = m x o (2-3) 

by its construction. Therefore by measuring the endpoint of Mt 2 distribution, one can 
measure the mass of the heavier of parent particles. Moreover, the LSP mass can be 
determined by measuring a kink position of the M™ 2 ax (m x ) as a function of the test mass 
m x [19,21,22]. Throughout this paper, we fix m x = m^o for simplicity. 

To determine masses of squark or gluino from the Mt 2 endpoint, decay products 
of sparticles should be correctly identified to i>j. Unless the event has relatively simple 

topology such as pp ^qq (jXi)O'Xi) or PP -> 99^ {j3Xi)(j3Xi), some of tlie j ets are 
mis-grouped into wrong V{. It has been proposed in [31,32] hemisphere algorithm is useful 
to define V{. It is defined as follows; 

1. Take two jets as the first seeds of the two hemisphere Hi and H 2 , Ji € Hi and 
J 2 £ H 2 : Ji is the highest pr jet, and J 2 is the jet (i) whose pn x AR(p(Ji),p(i)) is 
the largest in the event. 

2. Associate the other jets to the one of the hemispheres based on a distance measure 
d, so that jk belongs to Hi(i = 1, 2) if d(i, k) < d(j, k). Here d is defined as 

E (i) 

d(p£L,Pfc) = (4L - bhlil cosM hemi (2.4) 

l^hemi + ^) 

3. Take the sum of the jet momenta that belong to Hi and regard it as a new seed. 
Repeat the processes 2 and 3, while keeping the first seeds Ji and J 2 in the different 
hemispheres, till the assignment converges. 

After defining the hemisphere, we define the inclusive Mt2 using p Vi = ^2k e H.Pk- I n this 
paper we only use jets with pt >50GeV and |r/| <2.5. 

Alternatively, it was proposed to use MTGen which is the minimum of the Mt 2 s for all 
choices of jet combinations [18] 

M TGen = min M T2 (p Ha ,p H <x ,pf iss , m x ) (2.5) 

I<a<2™- 1 -1 

where a denotes a possible combination to split jets into two groups Hf and H^, and min- 
imization is taken over all possible 2 n_1 — 1 combinations. Furthermore, pH a = J2keH a Pk 
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and n is the number of visible objects. If all jets and leptons are decay products of A\ and 
A 2 , M T Gen < max{m Al ,m A2 }. 

At hadron collision, there are also particles coming from ISR in addition to those from 
sparticle decays. Transverse momentum of the ISR jet can be as large as those from squark 
and gluiino decays. The ISR jets in the visible systems v% and v-i smear the Mt2 endpoint 
significantly. 

In Ref. [26] a new definition of Mt2 that reduce ISR effect is defined, and mass de- 
termination based on the quantity is demonstrated for a pp — > ggX process followed by a 
gluino decay g — > jjxi, where X is a ISR partons. The ISR improved Mt2, Mr 2 (min), is 
defined as follows; 

1. Calculate Mr 2 («) which is a Mt2 calculated with a given grouping procedure but 
without involving the i-th. jet in pt order. In this paper, we take hemisphere algorithm 
after removing the i-th. jet for the grouping. 

2. Take the minimum of the Mr 2 (i) over the certain range of i, 

M T2 (min) = min M T2 (i). (2.6) 

i 

In Ref. [26], the minimization is taken up to the fifth jet and Mt2 is calculated by 4 jets 
because a gluino forced to decay into jjXi- The contamination above the endpoint is 
significantly reduced for Mr^min) distribution. In the paper, however, only the limited 
process was studied, and the techniques should be extended to the other SUSY processes, 
especially gluino-squark co-production because it is the dominant production process for a 
wide region of the parameter space (See Table [j].). 

A merit to use another definition of the Mt2 for gluino-squak co-production, M^ 2 (sub), 
has been discussed in Ref. [24]. It is defined as the Mt2 but the highest px jet is not 
included, namely, 

M T2 (sub) = M T2 (1). (2.7) 

In the parameter region where mg > nig, the subtraction of the highest pt jet tends to 
reduce the full-system of q-g into the subsytem of g-g or g-Xi (Xi denotes charginos or 
neutralinos.) effectively because squarks either decay into g or x% with high probability. In 
this case, we have 

M T2 (sub) < m-g. (2.8) 

It is shown that the endpoint of the iWr 2 (sub) correlates with gluino mass very well. By 
definition, M^ 2 (min) < Mj- 2 (sub), therefore the MT 2 (min) distribution may improve the 
sensitivity to the gluino mass by reducing the tail of the distribution. 

The parameter region with rriq < m g is not considered in Refs. [23,24,26]. For such 
a parameter region, gluino decays into jq, and the squark further decays into jxi- The 
jet from the squark two body decay has significant energy of the order of \pTj\ ~ niq/2 
for CMSSM like mass spectrum with nig ^> m^. The pt of the jet from the two body 
decay tends to be harder than that of ISR jet, therefore including the highest pt jet for the 
minimization for M7 i 2 (min) may remove a jet from the squark decay with high probability, 
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Table 2: Summary table of previous studies on inclusive Mti- 



making the distribution near the endpoint rather flat. Therefore, we introduce a variable 
M^ 2 od (min) that is defined as 

Mf 2 od (min) = mm M T2 (i) (2.9) 

for nig < nig. The quantities defined in Eqs. Q2.6| ), (E^) and ( |2.9| ) can be straightforwardly 
extended to the ones using MrGen variables. 

Table. ^ summarises the results of inclusive Mt 2 studies which have been done so far. 
The aim of this paper is to provide more comprehensive study of mass determination using 
the quantities defined above, in a wide range of CMSSM parameter space. 



3. Validation of hemisphere algorithm 

One of the difficulties to apply Mt 2 distirubiton for general SUSY processes is in defining 
appropriate two visible systems, v\ and v 2 . If v\ and v 2 are originated from decays A\ and 
A 2 , Ai — > ViXi, th- e Mt 2 (We call it Mr^true).) is bounded above by maxjm^ , niA 2 }- 
When there are several jets and leptons in an event, it is generally difficult to find such a 
"correct" assignment. 

Finding the correct assignment may not be necessary to reconstruct the endpoint. 
Instead, we may use a jet assignment algorithm that possesses the following properties: 

(i) Mt 2 calculated from the algorithm does not exceed max{m j 4 1 , niA 2 } significantly, 

(ii) A large number of events contribute to the endpoint region. 

Fig. H (la) shows the Mr2(true) distribution in parton level. As can be seen, the 
M^2(true) distribution has a sharp endpoint at max{niq,nig} = niq ~ 910 GeV. We gen- 
erate 60000 events at Point 5 with 7TeV proton centre of mass energy, and ISR is not 
included. We apply the following "minimal cuts" to reduce the standard model back- 
ground 1 

1. £^ iss > max(100GeV,0.2M cfr ), 

2. M cS > 500 GeV, 

o A rjets(parton) . . 
°" iv p T >100GeV - x ' 



1 Recent studies adopt much higher cuts to reduce backgrounds, but they do not reduce the signal 
significantly. 
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Figure 2: (la) — (Id): A/t2 distributions. (2a) — (2d): Distributions on a (Mt2(w^o) 



sign(— /8i2)vP^|) plane. (3a) — (3d): Distributions on a (M^fm^), sign(— /3^ 2 )-\/|/3j 2 |) 

plane. (4b) — (4d): Distributions on a (2iV wrong /A r part on, sign ( — j5f 2 ) Vw^iiT) plane. The cor- 
rect, random, Mtgcu and hemisphere assignments arc adopted in figures (ia), (ih), (ic) and (id), 
respectively, where i = 1 — 4. We generate 6 x 10 4 events at Point 5 with 7 TeV proton centre of 
mass energy. 



in addition, we also define 

4- CS>4 ("4 jet cut"), or A^<£Sg& > 2 ("2 jet cut"), 
for section || 
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Figure 3: A schematic picture of an assignment that provides a dangerously large Mt2 value in 
an event with \fi ^> Mi + M2. Dashed arrows represent the missing particles. Visible objects are 
grouped into X or Y. 

If Vi is not defined appropriately, Mt2 can exceed maxlm^, mA 2 ). Fig. |2| (lb) shows 
the Mt2 distribution where partons are randomly assigned into v\ or V2 so that each system 
has at least one constituents. We can see that the distribution has a large tail, and the 
endpoint structure is not seen. 

To see what kinds of events and assignments generate large Mt2 values, we show distri- 
butions with the random assignment on (Mt2, max(m«i,m„2)) and (Mt2, sign(— f3\ 2 ) v\M^\) 
planes in Fig.^ (2b) and (3b), respectively. Here, (3\ 2 = P r P 2 is an inner product between 
velocity vectors of A\ and A2 in the lab frame. From figure (2b), we can see that the Mt2 
is linearly dependent on the heavier of m v \ and m V 2- Thus, the assignments that provide 
large Mt2 also provide large m Vi . In addition, from figure (3b) we can see that the endpoint 
of the Mt2 distribution increases as increasing sign( 

-£12) VdW from to °- 8 - The tail of 
the Mt2 distribution is mainly caused by the events with large sign(— (3f 2 ) ylPvzl- 111 such 
events, A\ and A2 are highly boosted into back-to-back. 

Fig. ^ is a schematic picture of a wrong assignment that provides large Mt2- In the 
picture two initial particles are highly boosted into back-to-back, and daughter particles are 
incorrectly assigned into groups, X and Y. As can be seen, if two daughters, i and j, have 
different origin, the back-to-back boost makes their angle 9ij and momentum magnitudes 
large. This results a large invariant mass of the visible systems. 

For the events with large back-to-back boosts and two daughter particle momenta from 
the same sparticle makes small angle compared to ones from different origin. The idea of 
hemisphere algorithm is based on this observation, although the angle 6%j is modified into 
more elaborate distant measure, d. Interestingly, even if the algorithm fails to find the 
correct assignment, m Vi would not be too large because members in the same group have 
relatively small d. In the algorithm, the highest pr jet is assigned into a different group 
from the one contains the jet with the largest p^Ai?(p(Ji),p(z)). This also prevents from 
making m Vi too large. 

Fig. H (lc), (2c) and (3c) correspond to (lb), (2b) and (3b), respectively but hemisphere 
algorithm is used. From figures (ic) (i = 1, 2, 3), we can see that the tails in figures (ib) are 
removed significantly, and the endpoint structure at max{mA 1 ,niA 2 } is recovered. Fig. ^ 
(4b) and (4c) shows distributions on (2A r wrong /A'p arton , sign ( — j3\ 2 ) 1 j3'f 2 \ ) planes for the 
random and hemisphere algorithm, respectively. Here, N WIOQg is the number of wrong 
assignments. As expected, in events with large sign 

(-PlilVlfiwl ( ~ °' 4 )' the algorithm 
successfully selects the correct assignment with high probabilities. Compared to (4b), the 
improvement is significant. Note that large back-to-back boosts require 

> "Mi + mA 2 , 
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where y/I is centre of mass energy of colliding partons. An event sample used in Fig.|2] is 
at Point 5 with \/S = 7TeV, which has the largest sparticle mass scale and the smallest 
centre of mass energy of protons in our samples shown in Table. [|. The efficiencies of the 
algorithm would be better for the other samples, or at higher yS 

Another known algorithm is MxGen method. By its construction, this variable is 
smaller than or equal to Mr 2 (true) in event by event basis. Fig|2] (id) (i = 1 — 4) correspond 
to (ic) but MxGen method is adopted. From figure (4d), we can see that Mrcen can 
select the correct assignment in events with large sign ( — /3^ 2 ) y|^^| ( ?t 0.4) as hemisphere 
algorithm. This is because in such events all wrong assignments can provide larger m Vi 
than that from the correct assignment because of the back-to-back boost. 

Compared to hemisphere algorithm, MxGen distribution does not have a tail beyond 
maxjm^ , m,A 2 } in parton level without ISR. We fit distribution (lc) to two linear functions, 
/(M T2 ) = a±(M T2 - M* d 2 ge ) + b for sign(M T2 - My dge ) = ±1, and (Id) to a linear function, 
f(M TGen ) = c - dM TGen . The fitted endpoints is 912.2 ± 9.6 GeV for (lc). On the other 
hand, M^^ n = c/d are 906.8 ± 85.2 for (Id). The main source of the large error in 
the M^Qz n is poor statistics near the endpoint region. To improve the error we can fit 
the distribution away from the endpoint. However in this case, the central value become 
lower than the input squark mass. Both results are consistent with max{mq, m g } = rriq ~ 
910 GeV. 

As can be seen, the slope of the MxGen distribution is sharper than that of the Mt2 
distribution with hemisphere algorithm. This feature relies on the property that the MxGen 
can not exceed the mass of a sparticle. We will show that this property is not held for 
the MTG en (min) and the M™£ d n (min). In the rest of this paper, we adopt hemisphere 
algorithm unless it is explicitly stated, because both algorithms have similar performance 
but MxGen calculation is more computational intensive. In addition, as we will show later, 
hemisphere algorithm can be also applied to find out decay topology in the events. 

4. Parton level M^min) and M™ 2 od (min) distributions with and without 
ISR 

In section [2], we defined Mr 2 (min) and M™ 2 od (min) variables for gluino mass measurement. 
The idea of these variables is that the events above the Mr 2 (true) endpoint can be removed 
by subtracting the hardest ISR jet. Although those variables are designed for the events 
with a hard ISR jet, we first show parton level M^ 2 (min) and M^ od (min) distributions 
without ISR. SUSY processes are the mixture of the events with and without hard ISR 
jets. We therefore check that parton level M^nim) and M™ 2 od (min) distributions do not 
"collapse" near the endpoint by removing a parton from sparticle decays. 

Fig. U (left) shows Mr 2 (min) distributions, where red-solid and -dashed distributions 
correspond to Points 1 and 2 (nig > rriq), and blue-solid, -dashed and -dotted ones cor- 
respond to Points 3, 4 and 5 (rrig > m g ), respectively. Here we generate 60000 events at 
14TeV, and the minimal cuts defined in the previous section is applied. For Points 3, 4 
and 5, the endpoints of the distributions roughly agree with the gluino masses 636, 646 
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Figure 4: Left; 771^2 (min) distributions in parton level. Right; m™2 d ( mm ) distributions in parton 
level. Red-solid, -dashed, Blue-solid, -dashed and -dotted distributions correspond to Points 1, 2, 
3, 4 and 5, respectively. We generate 60000 events at 14TeV. 
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Figure 5: Schematic picture depicting decay patterns for various mass spectra. A, B, C and D 
correspond to uiq ^> nig, rriq > (m^ +m t ) > m§, rriq > wig > (mj +mt) and m g > rriq, respectively. 
Red and blue triangles represent jets from squark and gluino decays, respectively. Dashed arrows 
represent the decays are not main modes. 



and 651 GeV, respectively. On the other hand, for Points 1 and 2, the endpoints of the 
distributions are about 100 — 200 GeV smaller than the input gluino masses. 

The difference in the distributions at Points 1, 2 and those at Points 3, 4, 5 comes from 
the ordering on the gluino masses. Fig. ^ shows various mass spectra and corresponding 
decay patterns. Points 5 and 3, 4 correspond to types A and B in Fig. ||, respectively. A 
squark dominantly decays to jg, and a gluino decays to jjxi m the region. Although there 
is a large mass hierarchy between gluino and weak gauginos, px of jets from the gluino 
decay are relatively mild because of the three body decay. Consequently, the events tend 
to have multiple jets with modest px- 

On the other hand, mass spectra at Points 1 and 2 correspond to type D in Fig.|5[ A 
squark decays to two body final state jxi producing a high pt jet for the mass spectrum. 
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Figure 6: Left; parton level mT2(min) (black) and mxGen (min) (red) distributions at Point 5. 
Right; parton level m™2 (min) (black) and m™Qg n (min) (red) distributions at Point 1. 



Because a gluino decays to jq, SUSY events usually contain two high px jets from the 
squark decays. If we subtract one of the high pj- jets from the event, the remaining system 
loose too much energy and calculated Mt2 can not reach at the gluino mass. 

A variable M™ 2 od (min) is designed to reconstruct the gluino mass in type D spec- 
trum. In the definition, we do not subtract the two highest pt jets. Fig.^J (right) shows 
M™ 2 od (min) distributions in parton level, where the colour and line schemes are the same as 
the LHS figure. The endpoints for Points 1 and 2 are consistent with the input gluino masses 
612 and 629 GeV, respectively. On the other hand, for Points 3, 4 and 5, the endpoints 
are about 100 GeV larger than the gluino masses. The events above the input gluino mass 
are mostly g-q production events. This suggests that -M^mm) and M™ 2 od (min) should 
be used for appropriate mass spectrum and decay pattern. Namely we need to know the 
mass ordering of gluino and squark in the gluino measurement. We will discuss this point 
in the next section. 

M™^ d n (min) and M^ d n (min) variables can be defined analogously to M™ 2 od (min) and 
M™ 2 od (min), respectively. The left figure in Fig. |6| shows comparisons between MT2(min) 
(black) and MTc7 en (min) (red) distributions at Point 5. Both the distributions have end- 
points near the input gluino mass but also tails which mainly come from q-g production 
events. Similarly, the right figure in Fig. ^ shows M™ 2 od (min) (black) and M™^ d n (min) 
(red) distributions at Point 1. The endpoints are given by the input squark mass. At Point 
1 (rrig > rriq), g-g production events can produce M™ 2 od (min) and M™^ d n (min) larger than 
the squark mass. However the cross section of this production process is subdominant 
and the tails are tiny. To reconstruct the gluino mass, inclusion of ISR is necessary as 
we will show in the next section. We may slightly improve the distributions by adopt- 
ing M^° d (min) instead of M^™ od ^(min), although it would be computational intensive 
especially for the events with a large number of jets. 

We now discuss effect of hard initial state radiations associated with hard processes at 
LHC. It was explicitly shown in Ref. [26] that 

1. The hardest ISR parton in g and q production processes could be harder than partons 
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Point 1 


Point 3 


Point 5 


Point 3' 


a(ggj : ) (inclusive) /cr(gg) (exclusive) 


1.48 


1.61 


1.72 


1.42 


a (9Qj) (inclusive) /cr(gq) (exclusive) 


0.69 


0.77 


0.84 


0.69 



Table 3: Ratio between the number of events with and without a hard ISR parton. See text for 
the detail 

coming from g three body decay g — > jjxi- The ISR causes endpoint smearing of 
the inclusive Mt2 distributions. Significant events appears beyond the Mr2(true) 
endpoint for g-g production when the Mt2 is calculated from the four highest pr 
jets. 

2. The My2(min) distribution is less affected by ISR, because the leading ISR effect can 
be removed. It is explicitly shown that M™ 2 ax (min) ~ m g for the process pp —> gg 
when a gluino is forced to decay through three body decay mode g — > jjxi- 

The ratio of events with a hard ISR to those without any hard ISR is summarized 
in Table^J. Here we show the ratio of the number of events, iV(inclusive ggj or gqj 
sample)/A r (exclusive gg or gq sample) calculated by Madgraph/Madevent [55] with Pythia 
parton shower [56]. Here, the "exclusive sample" is generated from g-g or g-q process with 
parton shower, but if parton shower is resolved, namely if a cluster of the partons is isolated 
from the initial state with more than a certain /ct distance, the event is rejected. On the 
other hand the "inclusive sample" means the events generated from pp — > ggj or gqj matrix 
elements with full parton showers, where j denotes gluon or quark. Again, if the event does 
not have a resolved parton cluster, the event is rejected. Altogether, there is no overlap 
between exclusive and inclusive samples. The parton level distribution of the "matched 
sample" is therefore the sum of the events of the inclusive sample and exclusive sample. 
The kx cut off scale is chosen so that there are no discontinuity in the total distribution. 2 
As it has already discussed in Ref . [26] , the fraction of the inclusive events in ggj matched 
samples is always higher than that in gqj matched samples. This is because the difference 
in colour factors between g and q. 

In Fig. ^ (left) we show the px distribution of the ISR parton at Point 5 (a solid line). 
The pt is 200 GeV in average. The pt distributions of the ISR are roughly the same for all 
model points in this paper. On the other hand, the px distributions of the partons from 
squark/gluino decay are rather model parameter dependent. 

In the same figure we show the pt distributions of the highest pt partons from g decays 
for the ggj matched sample at Point 5 (thick dotted) and at Point 3' (thin dashed). The 
gluino dominantly decays through three body final state g —> jjxi a t Point 5. At Point 
3' the mass spectrum corresponds to type C in Fig. ^J, where rrig > (m^ + rat). In this 
type of spectrum, a gluino dominantly decays into two-body final state g — > tit, and the 
top and stop further decay to lighter particles. Thus, the highest px jet is relatively soft 

2 The kr cut depends on the shower algorithm. We use pythia pr ordered shower. The kr cut is 60 GeV 
and the pythia shower scale is 100 GeV. 
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Figure 7: Left: the pr distribution of ISR for ggj process at Point 5 (solid), the highest pt decay 
products at Point 5 (thick dotted) and Point 3' (thin dotted). Right: ratio of the pr of the third 
and first pr partons from gluino decays at Point 1 (solid) Point 5 (thick dashed) and Point 3 '(thin 
dashed) for the same process. The distributions are scaled from 20000 ggj matched events to 1000 
events at each point. 
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Figure 8: The px order of the ISR parton among all partons of the inclusive ggj sample. Left, 
center and right figures are at Points 1, 5 and 3' respectively. 



compared to that of Point 5 and the scale of px is close to that of ISR. On the other hand, 
the distribution of the highest px jet at Point 1 is similar to that at Point 5. 

Not only the distribution of the highest px jet, but also the distribution of the other 
jets could be different. In Fig.0 (right), we plot the distributions of ptz/pti at Points 
1 (solid), 5 (thick dashed) and 3' (thin dashed). At Point 1, the highest px jet takes 
significant part of the total energy of the event because q — > q\i modes dominate. While 
at Point 5, four partons from the gluino three body decays have roughly the same order of 
px- At Point 3', the tendency is even stronger. 

In Fig.||, we plot the distribution of zisr where zisr is the order of an ISR parton 
among all the px ordered partons in the ggj inclusive sample at Points 1 (left), 5 (middle) 
and 3' (right). Distributions for qgj inclusive sample are similar to those for ggj. 
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Figure 9: Parton level Mr2(niin)(a solid line) and A/™ 2 od (min)(a dashed line) distributions for the 
matched ggj at Point 1. The M™ 2 od (min) distribution has more events near the endpoint because 
of the events that ISR jets are correctly removed by the minimization. The distributions are based 
on 40000 matched g-g events generated by Madgraph. 



At Point 1, the average number of jets in the events are rather small, but about a 
half of the events has more than four partons with px > 50 GeV for g-g production. The 
probablity that the ISR parton becomes the 1st or the 2nd hardest parton is smaller than 
the one to be the 3rd hardest parton, because at this point a squark dominantly decays 
into Xi an d a high px quark. At Point 5, a gluino decays into three body final state jjx- 
The ISR parton is one among the five highest px partons in that case. Finally, at Point 3', 
the ISR tends to be the highest pt jet. This is because, a gluino dominantly decays into 
tt at this point. The average number of jets is large, and at the same time each of the jets 
from squark/gluino decay has relatively small px- 

We now consider the events that contribute to the endpoint of the Mr^min). For 
exclusive samples, the Mf2(min) endpoints should be smaller than those of the Mt2- Es- 
pecially at Point 1, removing one of the two highest pt partons from the system reduce the 
events near the endpoint too much as discussed previously. Instead of using the Mt2 (min) , 
we use M^ od (min) to obtain the endpoint. For the inclusive sample at Point 1, ISR tends 
to be the 3rd hardest parton of the event, and using the M^° d (mm) is justified also from 



M^2(min) for the ggj matched samples. The number of events near the input gluino mass 
612 GeV is quite small for the Mr2(min), while significant events remain for theM^ 2 od ( mm 
in the same region. Using the Mf° d (min) distribution clearly improves the sensitivity to 
the gluino mass. 

For Point 5, the minimization of MT2{i) over the five hardest jets is reasonable in 
removing the ISR jet, because the iisR distribution shown in Fig. ^ is flat for the five 
hardest jets. The endpoint of the M^2(min) distribution for inclusive ggj production 
would be near the input gluino mass because significant events has a hard ISR jet. The 
M^2(min) endpoint of the q-g exclusive sample would also be near the gluino mass. This 
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Figure 10: Parton level MyaCoiin) distributions of (top left) matched ggj sample and (bottom left) 
matched qgj sample calculated by Madgraph + Pythia. The corresponding jet level distributions 
obtained by HERWIG with detector simulation are placed in the right; (top right) g-g production 
and (bottom right) q-g production. 



is because if we remove the parton from g — > qj decay, the system becomes g-g, and the 
endpoint for such events is given as the input gluino mass. Therefore, by minimizing over 
the five hardest jets, a relation M™ 2 od,max (min) ~ m g is satisfied. The sum of the ggj and 
no-ISR qg productions provide the endpoint at nig. Additional ISR jets for g-q productions 
could lead the distribution that ends near the squark mass, which may smear the endpoint 
at nig. Fortunately, the fraction of gqj exclusive sample is less than a half of the total 
events at Point 5. 

Of course, the shape of the distributions near the endpoints can only be estimated 
by Monte Carlo simulation. In Fig. 10, we show the parton level M^mm) distributions 
of the matched sample (left) at Point 5 generated by Madgraph/Madevent with Pythia, 
and compare it jet level distributions generated by HERWIG with detector simulation. 
Here the jet level event is reconstructed by the AcerDET [57] with jet energy smearing 
of 5E/E = 50%/ '^J~(E(GeV)). Jet reconstruction of the AcerDET is replaced to the 
Cambridge-Aachen algorithm with AR = 0.4 using Fastjet [58]. 
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Figure 11: Flowchart of the ISR improved Mt2 analysis 



The top figures are the Mj^min) distributions for ggj production sample and the 
bottom figures are for gq production. At Point 5, the parton level endpoint for the g-q 
sample is around 700 GeV, which is consistent with the input gluino mass, 651 GeV and 
far below the input squark mass, 910 GeV. Note that the left and right distributions are 
roughly consistent although 1) Our Madgraph samples include matrix element contribution 
of the hardest ISR parton, 2) Our jet level samples obtained by HERWIG include all ISR 
effects in parton shower approximation. 3) HERWIG takes care spin correlation of matrix 
elements and all decay processes, while Pythia does not for sparticle decays. 



5. The results of jet level simulations 



In this section, we show jet level results and discuss mass parameter determination by 
M T2 (min) and M^ A (min). To make our discussion clear, we show the steps for mass 
determination in Fig. [ll[ We will see the importance of the proper choice of improved M?2 
depending on the observed event distribution. 



5.1 Mt2 distributions and effect of the ISR 

We compare the parton level Mt2 distributions with out ISR (left) to jet level ones 
(right) in Fig. 12, We generate the evens by HERWIG+AcerDET+Fastjet as discussed 
in previous secton. In the left figure, the distribution reduces rather quickly toward 
M mffl = max(mj, nig) and there are tails beyond M max due to mis-reconstruction of hemi- 
spheres. On the other hand, the jet level distributions at Points 1 and 3 do not show 
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Figure 12: Mt2 distributions at Points 1, 3 and 5 in parton level (left) and in jet level (right). 
We generate 60000 SUSY events at 14TeV. 



structures at neither the squark nor gluino masses. Especially, the number of events be- 
yond 750 GeV at Point 5 is not much different from that of Point 3 in jet level. 

In the previous studies [24], the endpoint of M?2 distributions are consistent with 
M max . The difference comes from the choice of the squark and gluino mass scale. In the 
paper, the gluino mass is around 800 GeV and the squark masses are taken from 881 GeV 
to 1561 GeV. As the squark and gluino masses increase, relative importance of the ISR 
jet reduces as the average px of the sparticle decay products increases linearly with the 
sparticle mass, while average px of the ISR jet does not increase so quickly. We take the 
scale of SUSY parameters rather light (nig < 700 GeV) in this paper so that squark and 
gluino masses are in the discovery region at the early stage of the LHC. 

The ISR effect reduces significantly when nig is increased for fixed nig. At Point 5, 
where the sum niq + nig is 350 GeV higher than that at Point 1, the endpoint is roughly at 
the input value, 910 GeV. The reduction of the ISR effect may be understood as follows. 
At this point, there are large mass difference between squark and gluino and the endpoint 
is saturated by My of the squark hemisphere, when hemisphere reconstruction is correct. 
The ISR jet has to be grouped into the "squark hemisphere" rather than the "gluino 
hemisphere" for the evnets contaminating beyond the endpoint, unless the pt of the jet is 
very high. 

We compare the parton level distribution and the jet level distribution at Point 5 
(niu L = 913 GeV) more quantitatively in Fig. [l^ (left). We found a linearly decreasing 
region of the Mx2 distribution is close to that of the parton level one. Fitting it by a 
linear function /(Mn) = a + bMx2 for both signal and tail regions, we find "kink" position 
at 887 GeV for jet level and 882 GeV for parton level with errors around 10% and 6% 
respectively, which roughly agree with the squark mass. 3 

In Fig.pl| (right), we also show the distributions at Point 4 (m^ = 837 GeV). The 
endpoint is more difficult to see, although slopes of the two distributions agree between 

3 Ax 2 / n -°- f ~ 1 for a11 fits - 
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Figure 13: Comparison between parton level Mt2 (solid lines) and jet level Mt-i (dashed lines) 
distributions at Points 5 (left) and 4 (right). We generate 60000 events and apply the 4 jet cut. In 
the left figure fitting lines are also shown. 

650 — 800 GeV. Extracting the mass scale from the distribution requires the detailed com- 
parison between Monte Carlo and data, and it is not the scope of this paper. Therefore we 
move on to the MT2(min) distributions which remove the hardest ISR jet efficiently. 

5.2 Events with Two high pt jets and proper choice of ISR improved M?2 f° r 
gluino mass determination 

In CMSSM like parameter region, event topologies are significantly different between the 
cases of rrig > rriq (Points 1, 2, I' and 2') and rriq > rrig (Points 3, 4, 5, 3', 4' and 5'). As we 
discussed in section |4|, in rriq > rrig region, squarks mainly decay to jg, and the j\i mode 
is subdominant. On the other hand, in rrig > rriq region, q — > j\i is the main mode and 
events contain two high pr jets. It is likely that the two hard jets are selected as the seeds 
of hemisphere reconstruction and thus assigned into different groups. Because each group 
has only one hard object, the masses of both visible systems are expected to be small. On 
the other hand, in rriq > rrig region, each visible system may have two modest jets from 
gluino three body decays, and the mass scale of visible systems represent the gluino mass. 

As discussed in the previous sections, the number of events beyond mm(rrig, rriq) will 
be significantly reduced for M^^inm) distribution. However, the shape of the distribution 
depends on the mass spectrum, and decay pattern. Especially the M^mm) distribution 
would be flat near the expected endpoint at Pointl I and 1', and we should use M™ 2 od (min) 
distribution to obtain the gluino mass. Another way to say, to obtain the correct gluino 
mass, we need criteria to use M™ 2 od (min) instead of Mr2(min) from experimental data. 

The fraction of events which survive under the 2 jet cut is the one of the key observation 



to make the choice. In Fig. [14|, we show the distribution of the events with at least two jets 
with pt > 200 GeV. The number of events after the cut is large at Point 1 (a black dotted 
line), compared to those at Points 3 — 5. As discussed already, this is because the q — > jxi 



mode dominates at Point I. If we see the excess, we should use M™ 2 od (min). Structure in 
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Figure 14: The jet level Mt2 distributions at Points 1 (black dotted), 3 (black solid), 4 (blue 
oslid) and 5 (red dashed) for the events with at least two jets with pt > 200 GeV. The sharp edge 
at Point 1 is mostly consisted by the events with 77,50 <4. The distribution at Point 2 is similar to 
that at Point 1 and is not shown. 



Mt2 distribution under the 2 jet cut might also be useful for determination of ordering of 
squark and gluino masses. We observe a sharp endpoint at the true squark mass (600 GeV) 
at Point 1. Although the number of the events after the 2 jet cuts is rather small, such 
structure also exists at Point 3. This is due to the squark decay into electroweak inos with 
significant branching ratio at this point. We will discuss about mixed use of M™ 2 od (min) 
and Mj^min) at Points 3 and 3' in the next subsection. 

Alternatively, one can use hemisphere mass m v to estimate fraction of the events that 
have gone though q — > j\ decays. Red (Blue) distributions in Fig. 15 (1) — (5) show 

the jet level m v distributions for Points 1 — 5 (1' — 5'), respectively. Here and 

(2) 

m\, are superposed in the distributions. The shape of the distributions are very different 
between m g > m g and m g > m g cases. In m g > rriq region (See Figs. |l~5| (1) and (2).), 
the distributions have a sharp peak at m v = 0, and the fraction of the other events are 
small. On the other hand, in rriq > m g (See Figs.|i~5| (3), (4) and (5).), the distributions 
have another peak around rn g /2, which is the contributions from g — > xjj mode. The 
fraction of the events with m v ~ becomes small. This suggests that the shape can give us 
information of the ordering of gluino and squark masses. We can also study the hemisphere 
mass after removing the jet z m i n , m^\i m - m ) and m^\i m - m ), where M^rnm) = M/^^mm)- 
The fractions of the events with mm(m^ (i m j n ), {i m in)) < 50 GeV are 73%, 45% and 
35% at Points 1, 3 and 5, respectively. 

Based on the discussions above, we decide to use Mf 2 od (min) for Points 1', 2', 1 and 
2, and Mr2(min) for 3' — 5' and 3 — 5 to determine the gluino masses. The jet level 
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Figure 15: Jet level m v distributions for the hemisphere reconstruction. The m%P 
superposed. Red (Blue) distributions in figures (1) — (5) correspond to Points 1 - 
respectively. 
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Figure 16: Left) M™ 2 od (min) (Point l':black) and Mr2(min) distributions (Points 3':red, 4':blue, 
5':purple). Right ) M™ 2 od (min) (Point 1: black ) and MT2(min) distributions (Points 3:red, 4:blue, 
5:purple). The mass of gluinos are from 522 — 558 GeV at Points 1' — 5' and 612 — 651 GeV at 
Points 1 — 5, respectively. 

distributions are shown in Fig. [0| The left figure shows the distributions at Points 1' — 
5' and the right figure show the distributions at Points 1 — 5'. The fomar have endpoints 
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point 


a 


5a 


b 


5b 


a/b 


gluino mass 


1 


2045 


84 


3.42 


0.16 


598 


612 


2 


2103 


107 


3.38 


0.20 


623 


620 


3 


2285 


64 


3.86 


0.13 


592 


635 


4 


1288 


57 


2.10 


0.10 


612 


646 


5 


1178 


82 


1.90 


0.13 


620 


651 


r 


1592 


70 


2.98 


0.16 


534 


522 


2' 


1571 


68 


2.95 


0.15 


533 


528 


3' 


1842 


67 


3.66 


0.15 


502 


541 


4' 


1765 


69 


3.37 


0.15 


524 


550 


5' 


1514 


81 


2.77 


0.17 


547 


558 



Table 4: The results of linear fits of the distributions in Fig. 16 near the endpoint. The fitting 
function is a simple linear function / = a—bMx2- For Points 1, 2, 1' and 2', M™ 2 od (rnin) distributions 
are used, and Mj^min) for the others. 



around 500 GeV while the latter is around 600 GeV. 

As we have seen in the parton level distribution, it is not good idea to focus on the 
events too close to the endpoint when we use the hemisphere algorithm to define Mti- 
Unlike Mxcen-, the algorithm tend to preserve the number of events near the endpoint, 
but there is a tail beyond the true endpoints because of the hemisphere misreconstruction. 
Therefore we fit the distributions to linear functions from the bin with the hight half of 
the maximum to the bin at the 1/10 ~ 1/6 of the maximum. The results of the fitting 
functions, f = a — bM^ od \min), are listed in TableQ. The errors are typically 7% and 
section of the line b/a are consistent with the corresponding gluino masses. The fitting 
function adopted here is too simple, leading xV n d.o.f are around 2. 

In the previous paper, we have proposed Mj^sub) to determine the gluino mass when 
rrig < rriq. The endpoint of M^Cnim) distribution agrees with the gluino mass better then 
My2(sub). The minimization among the possible assignment of jets from q —> gq decay 
among high px jets ensures MT2(min) distribution in qg production is below m g when 
there are no significant ISR jet. We show MT2(min) (red lines) and M^2(sub) (blue lines) 
the distributions at Points 5' (left) and 5 (right) in Fig. ^ with fitting curves near the 
endpoint. For comparison, parton level Mt% distributions for g-g samples without ISR jets 
are also shown by purple lines. The M^2(sub) distributions suffer higher tails and tend to 
overshoot the true endpoints. Exceptional case is Point 3', which we will discuss separately. 

5.3 Squark mass determination 

In this section we discuss squark mass determination from the experimental data. In 
Fig. 18, we show the distributions of M^ od (min) at Points 1' (left) and 1 (right) after 2 jet 



cut. Under the cut, the q-q production events are significant part of the survived events, 
because events with g — > it tend to be rejected. The structure due to the q-q production 
events can be seen in the total distribution. The endpoints of q-q parton level distributions 
(the thick dotted lines) are consistent with the positions of the edge structure in the jet level 
distributions (the solid lines), on the other hand g-g distribution (the thin dashed lines) 
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Figure 17: Mr2(min) distributions (red lines), Mxa(sub) distributions (blue line) and parton level 
Mt2 distributions for g-g production (purple lines). The left and right figures correspond to Points 
5' and 5, respectively. Fitting curves are also shown in the figures. 




Figure 18: Af™ 2 od (min) distributions after the two jet cut (solid lines), and parton level Mti 
distributions of qq (thick dashed lines ) and gg distributions (thin dashed lines). The left and right 
figures are for Points 1' and 1, respectively. 

are negligible. Taking M™ 2 od (min) reduces the contamination of q-g production beyond the 
edge. 

The jet distribution reduces very quickly near the endpoint. At Point 1', we fit the 
distribution by two linear functions. At point 1 we fit the distribution to a quadratic 
function bellow the kink and a linear function beyond the kink. The obtained kink position 
is at 483 (540.0) GeV at Point 1' (1). The values are very close to the input squark masses, 
and much smaller than gluino masses. Errors in the values are around 7%. 

Now we show distributions at Points 3' and 3. The squark masses is around 619 
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Figure 19: Mj^min) (black) and M™ 2 od (min) (green) distributions at Point 3' (left) and Point 
3 (right). The blue sold and blue dotted lines correspond to parton level gg and gq distributions, 
respectively. At Point 3 gluino decay into three body final state, while a gluino decays into a stop 
and a top at Point 3'. We also show a distribution of My2(sub) (red line) at Point 3' whose endpoint 
is slightly lower than that of M™ 2 od (min). 



(734) GeV at Point 3' (3). At these model points, the squarks are heavier than the gluino 
but the decay q — > jxi still has a significant branching ratio (See Table, [j]). Especially the 
fraction of the events survives under 2 jets cut at Point 3 is larger than those at Points 4 
and 5. Then there is a possibility to obtain the squark mass scale by fitting the M™ 2 od ( mm 
endpoint instead of Mr2(min). In Fig.^ we fit the events near the endpoint to a linear 
function. The endpoints of Mf?° d ( min) distributions are 557.4 GeV and 676.1 GeV with 
about 6 % of statistical errors at Points 3' and 3, respectively, higher than fitted endpoint 
of My2(min)(~ rrig) listed in Table||. Fits of the parton level Mt2 distirubtions of gq 
production process lead 601 GeV and 706 GeV at Points 3' and 3, respectively. 

Point 3' is a model point with special feature. At the point, gluino decays into tt 
dominantly. It is easy to see that g — > tt is the dominant gluino decay mode through the 
numbers of both jets and 6-jets in the events. The parton level Mt2 distributions at Point 
3' is suppressed near squark and gluino masses. This is due to the softness of average jet 
Pt- We only use the jets with pt > 50 GeV. On the other hand, the average pt of the 
jets are small so that some jets from gluino decays are not taken into account in the mass 
reconstruction. We may reduce the pt cuts but in that case the probability to assign jets 
into the wrong hemisphere becomes so large and the endpoints are smeared. Indeed, the 
endpoint of the parton and jet level distribution is systematically lower than the input 
squark mass. 

In the previous section, we have seen that the highest px jet is likely to be the ISR at 
Point 3' for g-g production. The main decay modes of squarks are q —> jg, and q — > jxi 
mode is subdominant. Therefore removing the highest pt jet might be equally useful to 
remove the hard ISR in qgj production. To check this, we also plot M^2(sub) distribution 
in the same plot (red line). Both M^2(sub) and M™ 2 od (min) have similar endpoints. 
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Figure 20: Mt% distributions under the 4 jet cut at Points 3' (red filled) and 5'(blue filled) and 
y/s = 7 TeV and those at 14 TeV( black solid line) . We use 30000 events at each model point. 



6. Expectation at y/s = 7 TeV and / Cdt = lfb -1 

Currently the LHC is operated at 7 TeV and will accumulate f Cdt ~ lfb -1 . Although 
the aim of this paper is to improve SUSY parameter determination in general hadron 
collider by using ISR improved inclusive Mt2, the event distribution at 7 TeV and lib -1 
of luminosity is also our interest. 



Again, we start with discussion of the Mt2 distribution at 7 TeV. Fig. 20 shows the 
Mt2 distributions at Points 3' (red filled) and 5' (blue filled). The figures are based on 
30000 events, but the expected number of events by the end of 7 TeV run is smaller (See 
Tabled). 

The smearing due to ISR is smaller than those at 14 TeV shown by the open histogram 
by black solid lines. The distributions beyond the true masses (shown by * marks and 
the lines) at 7 TeV are less than corresponding distributions at 14 TeV. Especially some 
structure remains around 600 GeV at Point 3'. At Point 5', Mt2 distribution is straight 
between 600 GeV to 800 GeV, and the endpoint obtained by fitting it to a linear function 
roughly consistent with the squark mass. The Mt2 distribution is correlated with squark 
mass and may be used for mass determination by fitting the distribution to the template 
based on Monte Carlo simulation data. 

In Fig. [2l|, we show the M|?° d ( min) distributions (bars) at Point 1' (Type D in Fig.|5[) 
with our 2 jet (left) and 4 jet (right) cuts together with the Mt2 distribution (green) under 
the same cut. In the figure we generate events corresponding to j C = lib -1 at 7 TeV. 
The coincidence of the endpoint with the input squark/gluino masses (478/522 GeV) are 
already visible. 

Fig. 22 shows the MT2(min) and M™ 2 od (min) distributions at Point 3', 3 and 5' which 



have mass spectraType A — C. At Points 3 and 5' production cross sections at 7 TeV 
are only around 1.7 pb and 2.5 pb and the distributions are under significant statistical 
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Figure 21: M T2 distributions at 7TeV and lfb -1 . Mf 2 od ( min) distirubtions with error bars at 
Point 1' for the 2 jet cuts (left), and the 4 jet cut (right). Green lines show Mti distribution under 
the same cut. 
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Figure 22: Mt2 distributions at 7TeV and lib -1 . The distributions of Mx2(min) (black bars) 
M™ 2 od (min) (red bars) and Mti (green lines) for Point 3' (left figure), 3 (central figure) and 5' 
(right figure). 

flactuation. 4 The Mt2 distribution ends around 800 GeV at Point 5, which reflects the 
input squark mass. Note that we have reason not to use M™ 2 od (niin) distribution for mass 
determination, because the number of the events which survives after 2 jet cut is smal. 
Therefore we do not show the M^ od ( min) distribution in the figure. 

7. Conclusion 

At the early stage of the LHC experiment, useful discovery channels are jets + -Ermiss 
channel and jets+ 1 lepton + Ermiss • The luminosity is rather low, so we want to measure 
sparticle nature from inclusive measurement rather than exclusive and clean modes. While 
Mt2 is useful kinematical variables in measuring parent SUSY partilce masses, an inclusive 
definition proposed in [23] is not protected from smearing due to ISR. 

4 For example, the distirbution at Point 5' ends around 600 GeV, but an endpoint structure at 550 GeV 
is found for 60000 signal events. 
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z-jets Mt2 


N(peak) 
of M hcmi 


IN(b-jets) 


determine d 


determine G 


Lj > b 


significant 
clear edge 


1 


few 


M^ 2 od (min) 
2-jets cut 


M|? 2 od (min) 
4-jets cut 


D > (j 

{g -> ft) 


mild 

edge 


2 


many 


Mf 2 od (min) 
2-jets cut 
or Mr2(sub) 


Mt2 (min) 
4-jets cut 


5 > G 
( 9 -> jjx) 


mild 

edge 


2 


few 


M^ 2 od (min) 
4-jets cut 


MT2(min) 


S<^G 


no clear 
edge 


2 


few 


M T2 


Mt2 (min) 
M T2 (sub) 



Table 5: Summary of the results obtained in this paper. G(S) denotes m§ and rriq respectively. 



In this paper we have proposed ISR improved inclusive Mt2 variables which might be 
useful to determine the squark and gluino masses separately. The modified M?2 variables 
discussed in this paper are Mj i 2(niin) and M™ 2 od (min). Endpoints of those distributions 
represent squark and gluino masses as summarized in Table ||. Steps to identify sparticle 
masses are shown in Fig. [ll]. The M™ 2 od (min) distribution is used for determination of 
gluino and squark masses for m g > nig case as the algorithm keeps the highest pt jets 
in the hemisphere which is likely comes from squark two body decays for rriq 3> m x case. 
Even if nig < m,, Br(q — > xq') is large enough so that M™ 2 od (min) is useful to determine 
nig. On the other hand, M^mm) is important to determine m g when m g < rriq and the 
decay of g is three body so that there are no dominant jets in the cascade decay. 

Errors in mass determination which we have obtained from the endpoint fits are typ- 
ically 6 — 10 % for 60000 generated SUSY events. On the other hand, we expect 4000 
(1000) events at 7TeV and / Cdt = lib -1 and 380000 (10000) events at 14TeV and 
J Cdt = 10fb _1 at Point 1 (5). Therefore, the model points with different gluino mass in 
Table|] can be clearly separated using the Mr2(min) or M^2° d (min) endpoints. We may 
identify an unique point in the squark and gluino mass parameter space up to the LSP 
mass uncertainty. This is very important meaning, because squark and gluino production 
crosss section is controlled by the masses up to small chargino and neutralino exchange 
contributions. We can cross check the predicted cross section to the observed number of 
events. In addition, the expected cross section would be significantly different if we change 
spin of the produced particles from scalar to fermion as is expected in Little Higgs models 
with T partity for example. Determination of the mass and decay pattern allows us to 
exclude/prove such models as well. 

There exist certain systematical errors through decay patterns of the SUSY parti- 
cles. In this paper, we study the mass spectrum that can be obtained in CMSSM. The 
procedure should be modified if mass hierarchies are significantly different from those in 
CMSSM. For example when wino or bino mass is much close to the squark mass, while the 
lightest neutralino remains much lighter than squark, we my expect moderate pt jets from 
neutralino/chargino cascade decays, even though rriq < m g , rather than prominent two 
high pt jets. Some of the squarks still directly decay into the lightest neutralino producing 
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a high pt jet if it is gaugino like. In such a case, although m g > nig, the distribution of 
the number of jets is close to that in rriq > m g case. Such classification of non-CMSSM 
models might be useful to reduce the systematics further. 
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